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ABSTRACT 



The solar system gas giant planets are oblate due to their rapid rotation. A measurement of the planet's 
projected oblateness would constrain the planet's rotational period. Planets that are synchronously 
rotating with their orbital revolution will be rotating too slowly to be significantly oblate; these include 
planets with orbital semi-major axes <; 0.2 AU (for Mp ~ Mj and M* ~ Jupiter-like planets 

in the range of orbital semi-major axis 0.1 AU to 0.2 AU will tidally evolve to synchronous rotation 
on a timescale similar to main sequence stars' lifetimes. In this case an oblateness detection will help 
constrain the planet's tidal Q value. 

The projected oblateness of a transiting extrasolar giant planet is measurable from a very high- 
photometric-precision transit light curve. For a sun-sized star and a Jupiter-sized planet the normalized 
flux difference in the transit ingress/egress light curve between a spherical and an oblate planet is a few 
to 15xl0 -5 for oblateness similar to Jupiter and Saturn respectively. The transit ingress and egress are 
asymmetric for an oblate planet with an orbital inclination different from 90° and a non-zero projected 
obliquity. A photometric precision of 10 -4 has been reached by HST observations of the known transiting 
extrasolar planet HD 209458 b. Kepler, a NASA discovery-class mission designed to detect transiting 
Earth-sized planets requires a photometric precision of 10 -5 and expects to find 30 to 40 transiting giant 
planets with orbital semi-major axes < 1 AU, about 20 of which will be at > 0.2 AU. Furthermore, part- 
per-million photometric precision (reached after averaging over several orbital periods) is expected from 
three other space telescopes to be launched within the next three years. Thus an oblateness measurement 
of a transiting giant planet is realistic in the near future. 

Subject headings: planetary systems — planets and satellites: general 



1. INTRODUCTION 

Extrasolar transiting planets are suitable for a variety 
of followup measurements while they transit their parent 
star. Many of the followup measurements would aim to 
detect additional planet components that block light from 
the parent star, such as moons, rings, or the planet's atmo- 
sphere. Recently, the atmosphere of the transiting planet 
HD 209458 b was detected in transmission in the neutral 
sodium resonance line (Charbonneau et al. 2002). A fol- 
lowup measurement of a transiting planet's albedo and 
phase curve would also be very useful in constraining at- 
mosphere properties since both the radius and orbital incli- 
nation (i) are known (Seager, Whitney, & Sasselov 2000). 
Here we describe an additional planet property that can be 
derived from followup measurements of transiting extraso- 
lar planets: the planet's projected oblateness. The oblate- 
ness is defined as (i? e — R p )/R p where R e is the planet's 
equatorial radius and R p the planet's polar radius. The 
solar system gas giants all have a significant oblateness, 
ranging from 0.02 to 0.1, due to their rapid rotation. Sat- 
urn's oblateness is so large it is obvious from whole-planet 
images. An oblateness detection would immediately tell 
us that a planet is a rapid rotator. 

A planet's oblateness signature on a transit light curve 
was first described in Hui & Seager (2002). The oblateness 
signature is expected to be small. However, recent obser- 



vations by HST/STIS (Brown et al. 2001) reached an un- 
precedented photometric precision of 10 -4 , and space tele- 
scopes MOST (Canadian Space Agency; launch date De- 
cember 2002), MONS (Danish Space Agency; launch date 
~ 2004), Corot (French Space Agency, European Space 
Agency, launch date 2004), and Kepler (NASA; launch 
date 2006) are designed to reach photometric precision of 
a fewxlO -6 by averaging measurements over several or- 
bital periods. Kepler, a wide-field transit search with the 
main goal of detecting Earth-sized planets orbiting solar 
type stars, expects to find 30-40 transiting extrasolar gi- 
ant planets. With the prospect of high photometric preci- 
sion, we investigate the transit signature of an oblate giant 
planet and the constraints on an oblate extrasolar planet's 
rotational period and tidal Q value. 



2. OBLATENESS AND ROTATIONAL PERIOD 

Here we review the relation between oblateness and ro- 
tational period (see e.g. Collins 1963; van Belle et al. 
2001), making the simplest assumptions about the struc- 
ture of the planets in question. The potential for a uni- 
formly rotating body has two terms, gravitational and ro- 
tational. The gravitational potential for a rotating, not 
too asphcrical, figure of equilibrium can be expanded us- 
ing the Legendre polynomials (see e.g. Danby 1962; Chan- 



drasekhar 1969), 
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where G is the Universal gravitational constant, Mp is the 
planet mass, R is the planet surface radius, Ji is a constant 
coefficient, and 9 is the angle measured from the planet's 
rotation axis. The gravitational moments, of which J 2 is 
an example, describe the structure and shape of a rotating 
fluid and are measured for the solar system planets by pre- 
cession rates of satellite and ring orbits and by spacecraft 
trajectories. The gravitational moments are not known for 
extrasolar planets, but they generally contribute small cor- 
rections to the gravitationalpotcntial. Here, we keep only 
the first term of equation (fy), the gravitational potential 
of a point mass. 

The rotational potential 4> r arises from the force of cen- 
tripetal acceleration. For a rigid body 1 rotating at angular 
velocity w about the z axis, 
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where the right hand side term is in spherical polar coor- 
dinates, and for axisymmetric rotation R = R(9). Hence 
the total potential <fi = 4> g (0) + 4> r {6) for a rigid rotator 
with a gravitational point mass potential is 
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We can use the fact that the total potential <fi is constant 
on the surface (Eddington 1926) to consider that (j> at the 
pole (8 — 0) and at the equator (0 — it/2) are the same. 
Then equation (0) can be written as, 
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where the polar radius R p corresponds to the rotation axis 
z, and R e is the equatorial radius. 

Then considering the definition of oblateness, rj = (i? e — 

Rp) / Rp; 



2GM F 



-V' 



(5) 



and the rotational period P = 2tt/lo is 
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(again, ignoring higher order terms in the gravitational 
potential). Considering a particle at r = R e and 8 = tt/2, 
GMpmj R 2 — moj 2 R e , and considering the relation be- 
tween the angular velocity uj and the oblateness r\ (equa- 
tion (|^)), the maximum possible oblateness is rj — 0.5, 
which corresponds to an ellipticity of e = 0.745. 

It is clear, considering equation (^), that a measurement 
of the oblateness 77, together with the measurement of Mp 
and Rp (here defined as Rp — yjR e R p ) from the transit 
photometry and radial velocity data, will yield the planet's 
rotational period, subject to two considerations. The first 
is that the (f> g is expected to deviate from a point mass 
potential, since the solar system gas giant planets have 
non-zero J2 and higher order moments. Secondly only the 
projected oblateness can be measured and hence only an 
upper limit on the rotational period can be obtained. To 
be precise, suppose the three-dimensional ellipticity is e, 
the projected ellipticity equals e(l — cos 2 a) where a is 
the angle between the axis of planet rotation and the line 
of sight to the planet (see Appendix A of Hui & Seager 
2002). To the extent the angle a is random, one expects 
on average that the projected ellipticity is half the actual 
value. The projected oblateness can be obtained from the 
projected ellipticity using the relation between r\ and e in 
equation (Q). Hereafter, whenever we refer to ellipticity 
e or oblateness rj, we mean its projected version, unless 
otherwise stated. 

3. TIMESCALE FOR TIDAL SYNCHRONIZATION OF 
ROTATION 

Planetary oblateness is caused by rapid rotation. Plan- 
ets that rotate synchronously with their orbital period will 
generally be rotating too slowly to be significantly oblate. 
Therefore, the tidal evolution timescale for synchronous 
rotation will limit the parameter space — mainly the orbital 
semi-major axis D — where oblate planets are expected to 
be found. Conversely, a detection of oblateness for a planet 
together with an estimation of the system's age could help 
to constrain the planet's tidal dissipation factor Q (we will 
refer to it as Qp here). The synchronous rotation timescale 
is (e.g. Goldreich & Soter 1966; Zahn 1977; Hubbard 1984; 
Guillot et al. 1996) 
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Note that with e as the eccentricity (hereafter called ellip- 
ticity to avoid confusion with eccentricity of the planet's 

3 2 , and so 



orbit) of an ellipse, R p = R e s/l 



primordial rotational angular velocity and u> rb its orbital 
angular velocity, R* and M* are the star's radius and 
mass, and G is the universal gravitational constant. Qp 
is inversely related to the tidal dissipation energy. For a 
Jupiter-like planet (with Jupiter's tidal Qp = 10 5 (Ioan- 
nou & Lindzen 1993) and Jupiter's current rotation rate 
u>i = ujupiter — 2-7r/9.92h) orbiting at 0.05 AU around a 
solar twin, t sync — 2 x 10 6 years. For HD 209458 b, with 
i?„ = 1.18i? , M* = l.O6M , R P = 1.42i?j, (Cody & 
Sasselov 2001; Brown et al. 2001), M P = 0.69Mj (Mazeh 
et al. 2000), t sync = 4 x Q P years. Because HD 209458 b 
has t S y nc "C t* (where i* is the age of the parent star), 
synchronous rotation is expected to have been reached. 

^^The gas giant planets have differential rotation. The slight deviation from rigid body rotation is due to convective redistribution of angular 
momentum. 
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The maximum oblateness can be estimated by consid- 
ering the centripetal acceleration for a given planet mass 
and radius that can prevent the body from flying apart, 
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system, the star's surface is represented by 



Note that planets of Jupiter's mass, Qp, and uii orbiting 
a sun- mass star with D <; 0.2 AU will be synchronous 
rotators. 

The tidal synchronization timescale has as its main un- 
certainty Qp and top It is interesting to compare t sync 
with the typical ages of stars (£*) with known extrasolar 
planets: one to a few billion years. Figure ||a shows t sync as 
a function of D and QpMp/M 2 to illustrate that planets 
with D ~ 0.15 to 0.4 AU have t sync similar to these typ- 
ical ages of stars with known planets. Of the parameters 
in equation (||), D, Mp, Rp can be measured for a transit- 
ing planet and Af* can be well estimated. Therefore Qp 
can be constrained with: an oblateness measurement as 
evidence for non-synchronous rotation; an assumption for 
uji] a planet with D such that t sync ~ t*. (Note that Qp 
can also be constrained by considering the tidal evolution 
time for tidal circularization, but in this case the values 
for Mp and Rp in specific cases are not known.) 

In addition to the planet's rotation rate, other orbital 
parameters will evolve as a consequence of tides raised 
on the planet by the star. The orbital circularization 
timescale (Goldreich & Soter 1966) and the timescale for 
co-planarity (defined as the coincidence of the planet or- 
bital plane and the stellar equatorial plane) are both longer 
than the synchronous rotation timescale (Rasio et al. 1966; 
Hut 1981). The planet's projected obliquity, (3 (defined 
here as the projected angle of inclination of the planet's 
rotation axis to the planet's orbit normal) is measurable 
for a transiting oblate planet (see §|]) that has an orbital 
inclination different from 90°. The planet's obliquity is 
also affected by tidal evolution. The timescale for tidal 
evolution to zero obliquity is the same as the tidal syn- 
chronous rotation timescale for a planet in a circular orbit 
(Peale 1999). For eccentric orbits the obliquity evolution is 
complex and coupled to other orbital evolution timescales 
(Peale 1999). However, a detection of obliquity together 
with other measured orbital and physical parameters for 
a planet in the tidal evolution regime will be useful to 
constrain the planet's evolutionary history. Note that a 
planet may escape evolving to synchronous rotation (or 
to a circular orbit and zero obliquity) in the presence of 
other nearby planets. HD 83443 (Mayor et al. 2002) is an 
interesting example described in Wu & Goldreich (2002). 

4. RESULTS AND DISCUSSION 

4.1. Transit Computation 

We want to determine the shape and magnitude of the 
transit light curve due to an oblate planet with different 
values of projected obliquity. Because the effect is small 
compared to a spherical planet transit, the transit light 
curve must be computed to high accuracy, <; 10~ 7 , in or- 
der to accurately calculate effects at the ~ 10~ 5 level. We 
use planetocentric coordinates and consider the projected 
oblate planet to be an ellipse. The planet's surface is de- 
scribed by the ellipse equation in polar coordinates, 
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where r is the distance from planet center and 8 is mea- 
sured from the projected orbit normal. In this coordinate 

2 Note that this definition of F is different from that in Hui & Seagcr 



R 2 = (r cos Q±D cos i cos u or b(t— t )) 2 + (r sinO+D sinuj or b(t— to)) 2 

(10) 

where u) or b is the orbital angular velocity, t is time, and 
to is the time of transit center. We consider the intersec- 
tion of the equations for the planet and star in order to 
determine (1) start and end of ingress and egress and (2) 
integration limits for the star area blocked by the planet. 

4.2. The Shape and Amplitude of the Oblateness 
Signature 



The oblateness signature on a transit light curve will 
be small. In order to quantify the shape and magnitude 
of the oblate planet's transit signature we compare an 
oblate planet transit light curve (T e m pse ) to a spherical 
planet transit light curve (T sp h ere ) for planets of the same 
projected area transiting the same sized star. We define 
F = \T sp here{t) — T e m pse (t)] where F is the flux difference 2 
between a spherical planet (projected to a disk) and an 
oblate planet (projected to an ellipse) transit light curve. 
The ellipse and projected sphere are always chosen to have 
the same total area i.e. the minimum of T sp h ere and the 
minimum of T e m pse are the same in the absence of limb 
darkening. Each transit light curve is normalized to the 
stellar flux. Figure |2|b shows the maximum value of F 
for different e and Rp/R*. The parameters i = 90° and 
j3 — 0° are chosen for this figure — as we will see, this gives 
the weakest oblateness signature, and so the contours in 
Figure ^b give conservative estimates of the required pho- 
tometric precision to detect oblateness. 

The transit light curve of a Jupiter-sized planet tran- 
siting a sun-sized star with orbital inclination i = 90° is 
shown in Figure ||a. The round-bottomed transit is due to 
limb darkening, in this case the solar limb-darkening value 
at 450 nm is used (Cox 2000). At this blue wavelength the 
sun is strongly limb darkened. The flat-bottomed tran- 
sit was computed neglecting limb darkening, or in a very 
non-limb-darkened color (e.g. at IR wavelengths, although 
solar limb darkening is minor at / band) chosen to illus- 
trate the range of possibilities. On the scale in Figure ||a, 
transits due to oblate and spherical planets are difficult 
to tell apart; we show their difference in subsequent pan- 
els. In Figure || and in this subsection we consider planets 
with the same area as Jupiter. Because Jupiter is oblate 
(with R e = 7.1492 x 10 7 m and R p = 6.7137 x 10 7 m), 
Jupiter's area corresponds to a sphere with effective ra- 
dius 6.9280 x 10 7 m. Figure Hb shows the flux differ- 
ence F = [T sphere (t) - T eUipse (t)} between the non- limb- 
darkened transit light curve of a spherical planet and an 
oblate planet with ellipticity e varying from 0.1 to 0.4. 
These transit light curves were computed for an orbital 
distance D = 0.2, but for the case of i = 90° the transit 
time scales as D 1 ! 2 . The flux difference F in Figure ||b is 
as high as 1.5xl0 -4 for an planet with ellipticity e = 0.4, 
corresponding to the oblateness of Saturn. 

The oblate planet transit signature can be explained as 
follows. At the start of ingress, when just a small fraction 
of the planet has crossed the stellar limb, an oblate planet 
(projected to an ellipse) covers less area on the stellar disk 
than a sphere (projected to a disk). Hence F (which is ac- 
tually the flux difference of the light deficit of a disk and 

(2002). 
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ellipse) is higher. As ingress progresses, the area of the 
star blocked by an oblate planet decreases compared to 
the sphere, until the ingress midpoint. At the midpoint of 
ingress, when the planet center coincides with the stellar 
limb, the area of the stellar disk that is blocked is roughly 
the same for an ellipse and a sphere (except for the devia- 
tion of the star limb from a straight line); hence the frac- 
tional difference at mid-ingress is close to zero. (F ~ at 
mid- ingress or mid-egress only for i = 90° and Rp <C i?*). 
The difference F also eventually reaches around towards 
the middle of the transit between ingress and egress. An- 
other way to consider the oblate transit signature is that 
for an oblate planet (with larger R e than a spherical one of 
the same area), the "first contact" will occur earlier, and 
the "second contact" will occur later, thus explaining the 
sign of the difference signal F. 

Figure ||c shows the same flux differences (Fs) as in Fig- 
ure 0b, but at the highly-limb-darkened color of 450 nm. 
In the presence of strong limb darkening, the maximum 
value of the oblateness signature at ingress/egress is lower 
than in the case of no limb darkening (Figure |b). This is 
because ingress and egress occur at the stellar limb where 
a given area takes out less luminosity than from closer to 
star center. Because of limb darkening the planet transit 
light curve is asymmetric around mid-ingress. While this 
is also true for a spherical planet, the effect is exaggerated 
for an ellipse. Limb darkening also causes a flux difference 
between an oblate and spherical planet when the planet is 
fully superimposed on the parent star; however this differ- 
ence is negligible compared to the oblateness signature at 
ingress and egress. 

In all of the above cases, Figure ||a-c, i = 90° and j3 = 
is assumed. An oblate planet with i < 90° and non-zero 
projected obliquity (3 has an interesting transit light curve. 
Figure |^d shows the transit light curve for different orbital 
inclinations (neglecting limb darkening). The dot-dashed 
line at i — 88.73° shows an almost grazing planet tran- 
sit. The transit light curves of an oblate planet and of a 
spherical planet are indistinguishable on the scale of Fig- 
ure ||d. Figure |e and Figuregjf show F, the flux difference 
in the transit light curve between a spherical planet and 
an oblate planet with j3 = 45° for various orbital inclina- 
tions. Note that ingress and egress last longer for transits 
with larger cosi, making the oblateness signature easier to 
measure. 

The ingress and egress oblateness signatures are com- 
pletely symmetric about mid-transit (i.e. F(t) = F(—t)) 
for projected obliquity — 0° (cases shown in Figure ||c). 
The oblateness signature will be more easily measured 
from a transit light curve that is asymmetric about mid- 
transit in its ingress and egress. Such an asymmetry will 
occur in the case of an oblate planet with a non-zero 
obliquity and with orbital inclination different from 90°. 
The asymmetry is apparent in comparing Figures ^b and 
|^e. Figures Si and b show the asymmetric flux differ- 
ence F asvm = [T e u ipse (t) - T eUipse (-t)] for different cases 
of i and (3. The oblate planet considered in Figure || has 
e = 0.4 and the same area as Jupiter. Reaching ampli- 
tudes of ~ 2 x 10~ 4 these are close to the current best 
precision photometry (from HST). The signatures in Fig- 
ure ^ can be considered maximum signatures using solar 
system values. 

In summary, there are two oblateness signatures one can 



look for: 1) a detailed difference in light curves between 
a transit by a sphere and a transit by an ellipse, around 
ingress and around egress; and 2) an asymmetry of the 
light curve between ingress and egress in the case of non- 
zero projected obliquity (3 and i ^ 90°. and e are not 
degenerate and can both be extracted from a fit to the 
light curve. When both signatures can be detected, one 
can obtain constraints on obliquity as well as oblateness. 

We expect to be able to distinguish the oblateness tran- 
sit signature from other similar magnitude effects during 
transits. Planetary moons could affect the planet tran- 
sit light curve but, due to their size, shape, and projected 
separation from the planet, would not mimic planet oblate- 
ness (although they could still cause a transit light curve 
asymmetry). Planetary rings will be projected into an el- 
lipse unless the inclination of the planet's ring plane to 
the line of sight is either 0° or 90°. With a few parts- 
per-million photometric precision, most cases of planetary 
rings should not be confused with the oblateness signa- 
ture because the rings have a relatively large radius and a 
slightly different transit signature than an oblate planet. 
A planet mass measurement (by radial velocity observa- 
tions and by the transit to give i) would likely alleviate 
any confusion due to transiting rings (which would make 
the planet look much larger) due to the giant planet mass- 
radius relationship (Guillot et al. 1996). Atmospheric 
lensing (Hui & Seager 2002) due to a spherical planet may 
be confused with the oblateness signature. However at 
wavelengths where the planet's atmosphere is strongly ab- 
sorbing atmospheric lensing is not present. In addition the 
transit asymmetry due to a non-zero (3 is not reproducible 
by atmospheric lensing. Confusion at the 10~ 5 level from 
other signals such as star spots, large planetary features, 
etc. are unlikely and are discussed in Hui & Seager (2002). 

4.3. HD 209458 b: The Only Known Transiting 
Extrasolar Planet 

Recently Brown et al. (2001) used HST/STIS to mea- 
sure the transit light curve of HD 209458 b, the only cur- 
rently known transiting extrasolar planet (Charbonneau 
et al. 2000; Henry et al. 2000). The resulting photomet- 
ric precision is at a level of 10~ 4 . This level of precision 
is not high enough to measure the expected oblateness of 
HD 209458 b. The close-in extrasolar giant planets like 
HD 209458 b are not expected to be significantly oblate 
due to slow rotational periods resulting from tidal locking 
to their orbital revolution. Assuming synchronous rotation 
an orbital period of 3.5 days corresponds to an oblateness 
of 0.0018 (e = 0.05), considering HD 209458 b to be a 
rigid rotator with M P = 0.69Afj (Mazeh et al. 2000) and 
Rp = 1.42 (the most recent derivation from Cody & Sas- 
selov 2001). At an orbital inclination of 86.1° the flux 
difference from a transit light curve of a sphere has a max- 
imum amplitude of only a few xl0~ 6 . If HD 209458 b 
has a non-zero projected obliquity (/3), the asymmetry of 
ingress to egress could be as high as 8xl0 -6 . Such an 
obliquity is not expected, however, due to tidal evolution. 
Obtaining an upper limit for e and (3 for HD 209458 b by 
a fit to a high-precision transit light curve would obser- 
vationally confirm the expected slow rotational period for 
HD 209458 b. 

4.4. Constraining the EGP Rotation Rate 
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The four solar system gas giant planets have significant 
oblateness, with values listed in Table 1 (note that oblate- 
ness here is the actual oblateness in three dimensions, i.e. 
not a projected value). With the exception of Jupiter, 
they all have substantial obliquity with values that are 
consistent with random (3. The solar system gas giant 
planets have measured rotational periods from 9.92 h to 
17.24 h (see Table 1). The rotational periods as calculated 
from equation (^|) using the measured oblateness values are 
within 20% or less of the actual rotational periods. The 
estimated rotational periods are consistently low, due to 
the neglect of the higher order gravitational moments in 
the gravitational potential (equation (El)). 

For extrasolar planet transit measurements only the pro- 
jected oblateness can be measured, and hence only an up- 
per limit on the rotational period can be obtained. With 
a measurement of the rotation rate together with the pro- 
jected oblateness (3 could be constrained. The rotational 
period could be directly measurable with the next gen- 
eration of radio telescopes for Jupiter-like planets with 
high synchrotron emission (like Jupiter) (Bastian, Dulk, 
& Leblanc 2000). Note that spectral line measurements 
may give information on winds in the planet atmosphere 
rather than the planet's rotation rate. 

5. SUMMARY 

An oblate planet and a spherical planet of Jupiter's 
area orbiting a sun-like star will have different ingress 
and egress transit light curves by as much as a few to 15 



xlO -5 . A planet's oblateness is most easily detectable in 
the case where the planet's orbital inclination is different 
from edge-on, and the planet's projected obliquity is sub- 
stantially different from zero. In this case the transit light 
curve is asymmetric between the ingress and egress (see 
Figure [l] and Figure ||) , a signature that can be detected 
by folding the light curve around mid-transit. The oblate- 
ness signature due to a planet with zero projected obliquity 
is symmetric in ingress and egress; the oblateness can then 
be measured using a model fit to the data (by exploiting 
the detailed differences in the light curves produced by a 
sphere and an ellipse; see Figure ||). The Kepler mission 
is expected to find 30 to 40 transiting giant planets with 
orbital semi-major axes < 1 AU, about 20 of which will 
be at D > 0.2 AU, and Kepler and three other space tele- 
scopes will reach part-per-million photometric precision. 
The future is promising for the oblateness measurement of 
extrasolar planets. 
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Planet 


Oblateness e 

R e — R p 
*- R, 


period (h) period (h) obliquity 
(measured) (calculated) 





Jupiter 


0.0648744 


0.34 


9.92 


8.23 


3.12 


Saturn 


0.0979624 


0.41 


10.66 


9.47 


26.73 


Uranus 


0.0229273 


0.21 


17.24 


13.83 


97.86 


Neptune 


0.0171 


0.18 


16.11 


14.07 


29.58 



Table 1 



Solar system giant planet oblateness, period, and obliquity from Cox (2000). Also listed is the period 
calculated from the oblateness by equation (g) (using mp and r e from cox (2000)). 




Fig. 1. — Definition of angles. The projected obliquity, /3, is the projected angle between the planet's rotation axis and the planet's orbit 
normal. The orbital inclination is i (where i = 90° corresponds to an edge-on orbit) and D is the semi-major axis. The planet in this schematic 
diagram has an oblateness equal to that of Saturn (10% which corresponds to e = 0.4; note that e in this paper is not orbital eccentricity but 
is rather related to oblateness, see §0.). For e > 0, i < 90°, and /3 > 0° the ingress and egress will be asymmetric, for example at points A 
and B. 
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Fig. 2. — Panel a: Synchronization timescales in years (contour lines) for different D, as a function of QpM P /M%. Typical ages of stars 
with known planets are a few billion years; Qp can be constrained from an oblateness detection with t sy nc of the same timescale. Note that 
a serious assumption in computing t sync is that loj = L0j up i ter = 2-7r/9.92 h. Panel b: The photometric precision (contour lines) required 
to detect a given ellipticity as a function of Rp/R„ (neglecting limtL darkening) . The parameters (3 = 0° and i = 90° were used; because 
cllipticity for /3 > 0° and i < 90° is much easier to detect (Figure S) the required photometric precision can be considered a conservative 
estimate. The corresponding oblateness, shown on the right y axis, is r\ = -2 — _ — 1. 
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-400 -200 200 400-400 -200 200 400 

Time (minutes) Time (minutes) 

Fig. 3. — A comparison of transit light curves for spherical and oblate Jupiter-sized planets orbiting a sun-sized star. Panel a: The 
normalized transit light curves (T) at orbital inclination i = 90° — on this scale a transit light curve from a spherical planet and an oblate 
planet are indistinguishable. The round-bottomed curve includes solar limb darkening at 450 nm, whereas the flat-bottomed curved is 
computed neglecting limb darkening. Panel b: F = [T sphere (t) — T e m pss (t)] (no limb-darkening) where T e u ip3e is for planets with projected 
ellipticity e = 0.1 (dotted line), e = 0.2 (short dashed line), e = 0.3 (solid line) and e = 0.4 (long-dashed line). Panel c: The same F as in 
panel b, but at the highly limb-darkened wavelength 450 nm. Panel d: The normalized transit light curve (T) neglecting limb darkening at 
orbital inclinations i = 90° (solid line), i = 89.47° (dashed line), i = 88.93° (dotted line), and i = 88.73° (dash-dot line). These values of i 
correspond to impact parameters (= cos(i)D/R t ) 0, 0.4, 0.8, and 0.95 respectively. Panel e: The flux difference F = [T aphere (t) — T e m pae (t)] 
at different i (line styles correspond to i in panel d) for a planet with e = 0.4 (that of Saturn) and /3 = 45°. Panel f: The same F as in 
panel e, but for solar limb darkening at 450 nm. The transit light curves and F in this figure were computed for D = 0.2 AU; for other 
orbital distances the time axis can be scaled by (.D/0.2AU) 1 / 2 (exactly for panels a-c and approximately for panels d-f). Note the asymmetry 
F(t) ^ F(—t) in the transit light curves in panels e— f (see Figure kj). 



8 



m 
i 

O 



X 



E 

in 
o 



30 
20 

10 



10 

20 
30 



E /X / \ / \ E 

L . ' x ■■ • \ \ J 


r h ; \ ! 

r Kj \i 


| \ 1 E 
' \ 1 -1 

\ / ! 

\y 1 



-400 



-300 -200 

Time (minutes) 



-100 







Fig. 



The asymmetric flux difference F a 



\T, 



ellipse 



(t) 



T, 



ellipse 



-t)] for a Jupiter-sized planet transiting a sun-sized star at 
D = 0.2 AU versus time from mid-transit (neglecting limb darkening). Panel a: F asyrn for e = 0.4 and i = 89.2° (corresponding to an impact 
parameter of 0.6), with /3 = 0° (solid line), f3 = 10° (dotted line), /3 = 20° (short dashed line), f3 = 30° (long dashed line), and /3 = 45° 
(dot-dash line). Panel b: Fasym for e = 0.4 and fi = 45°, with i = 90° (solid line), i = 89.73° (dotted line), i = 89.47° (short dashed line), 
and i = 89.20° (dot-dash line), and i = 88.73° (long dashed line) The i correspond to impact parameters (= cos(i)D/R„) 0, 0.2, 0.4, 0.6, and 
0.8, respectively. Note the solid lines in each panel show that there is no transit asymmetry for an oblate planet with /3 = 0° or i = 90°. 



